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We present a nonsingular bouncing cosmology using single scalar field matter with non-trivial 
potential and non-standard kinetic term. The potential sources a dynamical attractor solution with 
Ekpyrotic contraction which washes out small amplitude anisotropies. At high energy densities 
the field evolves into a ghost condensate, leading to a nonsingular bounce. Following the bounce 
there is a smooth transition to standard expanding radiation and matter dominated phases. Using 
linear cosmological perturbation theory we track each Fourier mode of the curvature fluctuation 
throughout the entire cosmic evolution. Using standard matching conditions for nonsingular bounc- 
ing cosmologies we verify that the spectral index does not change during the bounce. We show there 
is a controlled period of exponential growth of the fluctuation amplitude for the perturbations (but 
not for gravitational waves) around the bounce point which does not invalidate the perturbative 
treatment. This growth induces a natural suppression mechanism for the tensor to scalar ratio of 
fluctuations. Moreover, we study the generation of the primordial power spectrum of curvature 
fluctuations for various types of initial conditions. For the pure vacuum initial condition, on scales 
which exit the Hubble radius in the phase of Ekpyrotic contraction, the spectrum is deeply blue. For 
thermal particle initial condition, one possibility for generating a scale-invariant spectrum makes 
use of a special value of the background equation of state during the contracting Ekpyrotic phase. 
If the Ekpyrotic phase is preceded by a period of matter-dominated contraction, the primordial 
power spectrum is nearly scale-invariant on large scales (scales which exit the Hubble radius in the 
matter-dominated phase) but acquires a large blue tilt on small scales. Thus, our model provides a 
realization of the "matter bounce" scenario which is free of the anisotropy problem. 



PACS numbers: 98.80.Cq 

I. INTRODUCTION 

In this paper we construct a nonsingular bouncing cos- 
mology with an Ekpyrotic contracting phase. The model 
is free from the problems typically associated with 'mat- 
ter bounce' or pure Ekpyrotic models. Both the matter 
bounce (see e.g. pQ for a recent review) and the Ekpy- 
rotic scenario (see e.g [2] for an in-depth overview) were 
proposed as alternatives to inflationary cosmology as an 
explanation for the origin of the observed structure in the 
Universe. While these models are faced with serious chal- 
lenges we show that these challenges can be ameliorated 
by combining attributes of both scenarios. 

We begin with a discussion of the successes and prob- 
lems of the matter bounce model. Some time ago it was 
realized that curvature fluctuations which originate as 
quantum vacuum fluctuations on sub-Hubble scales and 
which exit the Hubble radius during a matter-dominated 
epoch of contraction acquire a scale-invariant spectrum 
[51 II]- In a number of toy models, it was shown that this 
spectrum often persists if the matter-dominated phase of 
contraction is continued in a nonsingular way to the ex- 
panding phase of Standard Big Bang cosmology. A non- 
singular bouncing model can be obtained by using Null 
Energy Condition (NEC) violating matter such as quin- 
tom matter [6] or Lee- Wick matter [Jj , by making use 
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of either a ghost condensate construction [5j, the ghost- 
free higher derivative gravity model of [9j[T0], or Galileon 
fields [HI [12] , and it arises in Horava-Lifshitz gravity [13] 
in the presence of spatial curvature [TU [TS] . A nonsin- 
gular bounce can also be obtained [TBHT5] in the context 
of mirage cosmology [HI [20] , by taking into account the 
effect of extra time- like dimension [2T], by making use 
of nonconventional K-essence model [221 123j , or in a uni- 
verse with nontrivial curved geometries [24H26] . In all of 
these cases it was found (under certain assumptions) that 
the spectrum of perturbations on scales relevant to cur- 
rent cosmological observations does not change its spec- 
tral index during the bouncing phase. Thus, such non- 
singular cosmologies provide an alternative to inflation 
for producing a scale-invariant spectrum of cosmological 
perturbations today. 

A difficult problem facing bouncing cosmologies is that 
the homogeneous and isotropic background cosmological 
solution is unstable to the development of radiation [27J 
and anisotropic stress. The latter instability leads to 
the famous BKL [2S] mixmaster cosmology, opposed to 
a homogeneous and isotropic bounce. 

In the Ekpyrotic scenario [25], the contracting branch 
solution is a local attractor [3U], similar to the accel- 
erating expanding cosmological solutions in inflation- 
ary models. In Ekpyrotic cosmology there is the as- 
sumed existence of a matter fluid with an equation of 
state w = p/p 3> 1 (p and p being pressure and en- 
ergy densities, respectively). The energy density in this 
fluid blueshifts faster than the contribution of anisotropic 
stress in the effective energy density. Thus, the fluid 
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comes to dominate during the contracting phase and pre- 
vents the development of the BKL instability. 

Within the context of a pure Ekpyrotic model, the cur- 
vature fluctuations on super-Hubble scales in the con- 
tracting phase are not scale-invariant; although, the fluc- 
tuation in the Bardeen potential [31j [j] which describes 
the metric inhomogeneities in the longitudinal gauge are 
scale-invariant (see e.g. [34-37 for discussions of this 
issue). The initial Ekpyrotic scenario [35] had a singu- 
larity in the effective held theory at the transition point 
between contraction and expansion. The question of how 
fluctuations transfer from the contracting to the expand- 
ing phase is non-trivial in this context (see e.g. [38] 
for a detailed discussion). There are prescriptions ac- 
cording to which the spectrum of fluctuations after the 
bounce is scale-invariant [39) . Also, there are typically 
entropy modes present in Ekpyrotic models. These en- 
tropy modes can acquire a scale-invariant spectrum by 
the same mechanism that $ acquires such a spectrum 
[10H33], and will then induce a scale-invariant spectrum 
for the curvature fluctuations by the usual mechanism 
of isocurvature modes seeding an adiabatic fluctuation. 
This leads to the so-called "New Ekpyrotic" scenario [H] . 
Scale-invariant entropy modes arise from fluctuations of 
the extra metric fluctuation modes in higher dimensions 
[45j , the setting in which the Ekpyrotic scenario was ini- 
tially proposed. However, in the context of the New 
Ekpyrotic scenario, one has to worry about instabilities 
during the bounce phase, an issue recently raised in |46j . 

The model we propose in this paper is free from the 
complications associated with entropy modes and the 
transfer of fluctuations through singular bounces. In this 
scenario, it is the curvature fluctuations which are scale- 
invariant. This spectrum is inherited from the phase of 
matter contraction which preceded the phase of Ekpy- 
rotic contraction. There are no low mass entropy modes 
leading to the complicated issues existing in pure Ekpy- 
rotic models. The cosmological bounce is nonsingular 
and hence the transfer of fluctuations from the initial 
contracting phase to the expanding period can be treated 
exactly. We find that the spectral index of the fluctua- 
tions on infrared scales relevant to current observations 
does not change during the bounce. The amplitude of the 
curvature fluctuations are boosted by a significant, scale- 
independent, factor JF, (for long wavelength modes). We 
demonstrate that our model is free from the instability 
problems raised in [41)] . 

In this model the background dynamics before and 
during the bounce are determined by a single scalar field 
(j) with non-trivial kinetic action of the form used in ghost 
condensate [37] and Galileon models [35] (see also [49] ) . 
We add to this field a negative exponential potential sim- 
ilar to what is used in the Ekpyrotic scenario. This po- 
tential leads to a phase of Ekpyrotic contraction. At 



high densities, the coefficient of the part of the kinetic 
term which is standard becomes negative, and this leads 
to a nonsingular cosmological bounce similar to what is 
achieved in the New Ekpyrotic scenario [42, 44J and in [5J. 
After the bounce, the coefficient of the kinetic term be- 
comes positive again, a period of kinetic-driven expansion 
sets in, lasting until the radiation and matter dominated 
periods of Standard Cosmology. 

We can imagine that the initial contracting phase of 
the universe mimics our currently observed universe, 
namely a state filled with regular matter and radiation. 
Since the kinetic energy density of <p increases faster 
in the contracting phase than both that of matter and 
radiation, eventually <fi will begin to dominate. Given 
this setup, our model provides a realization of the mat- 
ter bounce scenario of [TJ [3] H] which is free from the 
anisotropy problem which generically affects bouncing 
cosmologies. 

The outline of this paper is as follows: In the next 
section we introduce the model and discuss the cosmo- 
logical background dynamics, both analytically and nu- 
merically. In Section III we show how cosmological fluc- 
tuations evolve through the bouncing phase. The final 
section is reserved for discussion. 

We defined the reduce Planck mass by M p = 1/VSwG 
where G is Newton's gravitational constant. The sign of 
the metric is taken to be (+, — , — , — ). Note that we take 
the value of the scale factor at the bounce point to be 
as = 1 throughout the paper. 

II. MODEL AND BACKGROUND COSMOLOGY 

It is well known that, in order to realize a success- 
ful nonsingular homogeneous and isotropic bounce in 
a spatially flat (k = 0) Friedmann-Robertson-Walker- 
Lemaitre (FRW) universe in the framework of standard 
Einstein gravity, the Null Energy Condition for matter 
must be violated. This is because at the bounce point 
the time derivative of the Hubble parameter is greater 
than zero while the Hubble parameter itself is zero. This 
implies that the total energy density vanishes while the 
pressure is negative, i.e., that the background equation 
of state is w < —I. One way to achieving such a sce- 
nario is via a ghost condensate field cf> [SJ [35] [33] in which 
the Lorentz symmetry is broken spontaneously in the 
infrared, with the kinetic term for </> taking on a non- 
vanishing expectation value. However, this type of model 
suffers a gradient instability since the square of the sound 
speed of its perturbations becomes negative in the phase 
of ghost condensation Another mechanism of achiev- 
ing Null Energy Condition violation is to make use of a 
Galileon type field. The key feature of this type of field 
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is that it contains higher order derivative terms in the 
Lagrangian while the equation of motion remains second 
order, and thus does not lead to the appearance of ghost 
modes. Galileon fields have been used to construct emer- 
gent universe scenarios in which the universe begins with 
a quasi-static phase [50]. There have also been recent 
studies of Galileon models which yield cosmologies with 
a nonsingular bounce, but in which a space-time singular- 
ity of big rip type arises after the bouncing phase [TTHT2"] . 
In this section, we present a nonsingular bouncing solu- 
tion in terms of a single field having the desirable features 
of both the ghost condensate and Galileon-inspired mod- 
els. 



A. The Model 

The most general form of single scalar field La- 
grangian giving rise to second-order field equations in 
four-dimensional spacetime can be expressed as [511 152] 

C = K(<f>,X) + G(<j>,X)n<f) + L i + L 5 , (1) 

where K and G are functions of a dimensionless scalar 
field an d its canonical kinetic term 

X = . (2) 

The standard kinetic Lagrangian corresponds to K = 
X and all other terms vanishing. A more general form 
is K((f>,X) = A{$)B(X) with ,4(0) > used in "K- 
essence" models [55] . 

The other kinetic terms of include the operator 

□0 = <7^V„V^. (3) 

The terms L4 and L5 are higher order operators which 
are usually suppressed at low energy scales. Thus, in 
the present paper we will ignore them and focus on the 
first two terms adopting the form in the "Kinetic Gravity 
Braiding" (KGB) model of [55 (see also [55]). 

Variation of the above matter action minimally cou- 
pled to Einstein gravity leads to the modified Einstein 
equation: 

Tfi V = Mp (^R pv — — 9ij,u\ 

= (-K + 2XG^ + GxV.XV^)^ 
+(K, X + G >x n<t> - 2G )V M 0V„0 
-GxCV^XVW. + V^XV^) . (4) 

In the above formalism, and t x denote derivatives with 
respect to and X, respectively. 

For the model under consideration we choose: 

#(0, X) = M 2 V [1 - ,g(0)] X + PX 2 - V(<P) , (5) 

where we introduce a positive-definite parameter ft so 
that the kinetic term is bounded from below at high en- 
ergy scales. For g > 1 a ghost condensate ground state 



with X can arise. Note that the first term of K 
involves M£ since in the present paper we adopt the con- 
vention that the scalar field is dimensionless. We have 
also introduced a non-trivial potential V for 0. This po- 
tential is chosen such that Ekpyrotic contraction is pos- 
sible. In the specific model which we will discuss in the 
following, the scalar field evolves monotonically from a 
negative large value to a positive large value. The func- 
tion g{(j>) is chosen such that a phase of ghost condensa- 
tion only occurs during a short time when approaches 
= 0. This requires the dimensionless function g to be 
smaller than unity when \<j)\ 3> 1 but larger than unity 
when approaches the origin. 

The term G is a Galileon type operator which is con- 
sistent with the fact that the Lagrangian contains higher 
order derivative terms in 0, but the equation of mo- 
tion remains a second order differential equation. Phe- 
nomenologically, there are few requirements on the ex- 
plicit form of G. We introduce this operator since we 
expect that it can be used to stabilize the gradient term 
of cosmological perturbations, which requires that the 
sound speed parameter behaves smoothly and is positive- 
definite throughout most of the background evolution. 
For simplicity, we will choose G to be a simple function 
of only X: 

G(X) = -iX , (6) 
where 7 is a positive-definite number. 

B. The Cosmological Background Evolution 

We now study the equations of motion which result 
from applying the Lagrangian to the flat FRW uni- 
verse whose metric is given by 

ds 2 = dt 2 - a 2 {t)dx 2 , (7) 

where t is cosmic time, x are the comoving spatial coordi- 
nates and a(t) is the scale factor. Thus, the kinetic terms 
of a homogeneous and isotropic scalar field background 
become 

X = 4> 2 /2 , 
□0 = + 3i70 . (8) 

The dot denotes the derivative with respect to cosmic 
time t. The Hubble parameter is H = a/a. 
For this background the energy density is 

P = \m 2 p {1 - 5 )0 2 + \p4> A + 3 7 ^0 3 + V{4>) , (9) 

and the pressure is 

P = \m 2 p (1 - g)4> 2 + J/30 4 - 7 2 - F(0) , (10) 

by computing the diagonal components of the stress- 
energy tensor Q. 
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To derive the equation of motion for <j>, one can either 
vary the Lagrangian with respect to <f> or require that the 
covariant derivative of the stress-energy tensor vanish. 
This yields 



Vdp + v<t> + y = o 



(ii) 



where we have introduced 
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V = (l- g)M* + 6 7 #0 + m 2 + ^4> 



(12) 



V = 3(1 - g)M 2 p H + (9 7 # 2 - ^M 2 p g A )4> + 3pH<j> 2 



-^-g)i4> z 



9j 2 H<j) 4 3/37c/» 5 



2M| 



2M2 



(13) 



The coefficient V determines the positivity of the kinetic 
term of the scalar field and thus can be used to judge 
whether the model contains a ghost or not at the per- 
turbative level. The coefficient T> is an effective damping 
term. By keeping the first terms of the expressions of 
V and T> and setting g = one can recover the stan- 
dard Klein-Gordon equation in the FRW background. 
Neglecting the other terms is a good approximation when 
the velocity of <fi is sub-Planckian. 



C. The Nonsingular Bounce Solution 

We now consider a nonsingular cosmological bouncing 
solution. Here, we focus on homogeneous solutions. It is 
well known that the homogeneous trajectory of a scalar 
field can be an attractor solution when its potential is an 
exponential function. One example is inflationary expan- 
sion of the universe in a positive- valued exponential po- 
tential, and the other one is the Ekpyrotic model in which 
the homogeneous field trajectory for a negative exponen- 
tial potential is an attractor in a contracting universe. 
For a phase of Ekpyrotic contraction, we take the form 
of the potential to be 



V{4>) = - 



2Vn 



(14) 



where Vq is a positive constant with dimension of 
(mass) 4 . Thus the potential is always negative and 
asymptotically approaches zero when \(f>\ 1. Ignor- 
ing the second term of the denominator, this potential 
reduces to the form used in the Ekpyrotic scenario. 

For initial conditions we assume that <f> begins at an 
asymptotically large negative value. The force due to 
the potential induces motion towards the right, i.e. with 
<f> > 0. The equation of state of <f> matter is determined 
by the parameter q. For sufficiently small values of q 
(as will be shown shortly, the transition value is q — 
1/3), the equation of state is such that in a contracting 
universe the energy density in (f> increases faster than that 
in matter, radiation and anisotropic stress. Thus, the 
homogeneous trajectory is an attractor. As cf> approaches 
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FIG. 1. A sketch of the coefficient g and of the potential V as 
functions of the scalar field cf>. Depending on the value of g, 
the field space is separated into three regimes: the contracting 
phase with nearly canonical kinetic term (large negative field 
values) , the ghost condensate phase when g > 1 which occurs 
at field values close to 0, and the expanding phase with canon- 
ical kinetic term (large positive field values) , respectively. We 
find that the phase of ghost condensation leads to a nonsin- 
gular cosmological bounce. 



zero, the second term in the denominator of the potential 
becomes important and the field evolution departs from 
the Ekpyrotic trajectory. Since the potential is bounded 
from below the model has a stable vacuum state. 

To obtain a nonsingular bounce, we must make an ex- 
plicit choice of g as a function of <j). As we have discussed 
in the previous subsection, we want g to be negligible 
when \<j>\ 1. In order to obtain a violation of the Null 
Energy Condition after the termination of the Ekpyrotic 
contracting phase, g must become the dominant coeffi- 
cient in the quadratic kinetic term when <j> approaches 0. 
Thus, we suggest its form to be 



25o 



(15) 



+ e 



where go is a positive constant defined as the value of g 
at the moment of <fi — 0, which is required to be larger 
than 1 as will be discussed later. 

At this point we have fully determined the model. The 
potential V{4>) and the function g((p) are sketched as a 
function of (f> in Figure [T] The horizontal axis is the 
field value 4>, the vertical axis shows the values of the 
functions V (blue curve, negative-definite) and g (red 
curve, positive-definite). This figure is helpful in gaining 
a semi-analytic understanding of the evolution, the topic 
we turn to in the following subsection. 



D. Semi- Analytic Analysis 

As initial conditions for t <C —1 (in Planck units 
or units of the typical mass scale in the Lagrangian, 
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whichever yields a larger time), we assume a nearly ho- 
mogeneous universe filled with pressureless matter, radia- 
tion, and a homogeneous <f> field condensate with <f> <C — 1 
which is in a contracting phase. From Fig. [T] we see that 
in this phase g — > and (since the gradient of the po- 
tential also approaches zero) 4> <C M v . In this case the 



Lagrangian for 
form: 



approaches the conventional canonical 



M 2 

c^-g-d M m-v{<j>) , 



(16) 



This Lagrangian admits the Ekpyrotic attractor solution 
for <j) in a contracting universe: 



In 



q(l - 3q)M 2 



which yields an effective equation of state 



w c 



2 

1+ 3, 



(17) 



(18) 



where the subscript " c " denotes the contracting phase. 
Consequently, long before the bounce, the energy den- 
sity of the scalar field <fi evolves as that of a perfect fluid 
with a constant equation of state w c . Provided w c is 
the largest among the equations of state of all the mat- 
ter components in the universe, the contribution of 4> to 
the total energy density will become dominant, as is well- 
known in the Ekpyrotic model. If q < 1/3 we have w c > 1 
and in this case the importance of 4> increases also rel- 
ative to that of anisotropic stresses, demonstrating that 
the model is free from the BKL instability which plagues 
most bouncing cosmologies. 

As (f> accelerates towards = 0, the value of g will 
increase. If go > 1 (which we require), then at some point 
in time g will start to exceed the critical value g = 1 
and thus the sign of quadratic kinetic term in ([!]) will 
become negative. At that point, <f> will become a ghost 
condensate. The critical value of g which signals the 
onset of the ghost condensate phase is 



ff(0* 



1 



(19) 



There are two solutions of this equation which are 



ln2g 
P 

In 2g Q 
b g p 



and 



(20) 



The ghost condensate phase occurs for < <f> < (/)*+. 
During this phase the Null Energy Condition is vio- 
lated This allows for the existence of a nonsingular 
bounce. However, a nonsingular bounce also requires 



3 It is not sufficient to have p = which could occur for /3 
before ghost condensation sets in. 



that the energy density vanishes at the bounce point, 
which implies the following relation 



1 



M 2 V (1 - g B ) 



3 

4 ; 



V B =0 



(21) 



at the time t b when the bounce occurs (subscripts B refer 
to the bounce point). In our specific example, 4>b — 
and thus gs — go and Vb — Vq. Therefore, we find that 



. (g - l)M p 2 



3/3 




\2pVo 



(go - 1) 2 M4 



, (22) 



at the time of the bounce. From this result it follows 
immediately that a nonsingular bounce can only occur 
when go > 1 (otherwise there exists no real solution to 
Eq. (22)). 

Since it has a large positive velocity at the bounce 
point, 4> continues to increase after the bounce. Within a 
short time it will cross the second boundary of the ghost 
condensation region 0*+. At that point, the Lagrangian 
of the model recovers the canonical form and the universe 
enters a kinetic-driven phase of expansion. Note that al- 
though the potential approaches an exponential form, the 
scalar field does not approach the solution analogous to 
(17) which would be an attractor for positive values in 



a contracting universe. In an expanding space-time, this 
solution is a repeller rather than an attractor. Instead, 
the scalar field experiences a fast rolling phase with an 
effective equation of state the same as that of a stiff fluid: 



1 



(23) 



where the subscript " e " denotes the expanding phase. 
As a consequence, the energy density of the scalar field 
will dilute relative to that of conventional radiation and 
matter and then the universe will be able to enter the 
phases of the usual thermal history of Standard Big Bang 
cosmology. Following our qualitative analytical analysis 
we now turn to a numerical study of the background 
cosmology. 



E. Numerical Analysis 

Above, we have presented a heuristic discussion of 
how the Null Energy Condition can be violated in this 
setup. The key issue for the numerical analysis is to ana- 
lyze whether the Null Energy Condition can be violated 
smoothly and without any pathologies. In our model, we 
make use of the idea of ghost condensation to violate the 
Null Energy Condition. It is well known that any sin- 
gle field described by a K-essence type Lagrangian [53] 
will not cross the cosmological constant boundary (see 
e.g. Appendix of |56j). For a more general ghost conden- 
sate Lagrangian, violation of the Null Energy Condition 
is possible, but perturbations could become very large 
and force the background trajectory back away from the 
cosmological constant boundary. In addition, there are 
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FIG. 2. Numerical plot of the Hubble parameter H (verti- 
cal axis) as a function of cosmic time (horizontal axis). The 
main plot shows that a nonsingular bounce occurs, and that 
the time scale of the bounce is short (it is a "fast bounce" 
model). The inner insert shows a blowup of the smooth Hub- 
ble parameter H during the bounce phase. The background 
parameters were chosen as in Eq. (24 1. All numerical values 



are in Planck units M p . The initial conditions were chosen as 
described in the following figure. 



gradient instabilities [37]. To cure these, we take into 
account the effect of the Galileon operator G(X) = jX. 

In order to prove that the model is well-behaved 
throughout the entire cosmological evolution, we need 
to study the evolution of perturbations about the back- 
ground solution. We perform the analysis of cosmological 
perturbation in the next section. Here we will focus on 
the pure background dynamics. 

To illustrate that a nonsingular bounce can be achieved 
in our model, we numerically evolve the Einstein acceler- 
ation equation coupled to the field equation for <f>, impos- 
ing the Hamiltonian constraint equation to set the intial 
conditions. In the numerical computations we work in 
units of the Planck mass M p for all parameters. Specifi- 
cally, these parameters are chosen to be: 



V = 10- 



9o 



1.1 
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0.5 , p= 0.01 , q 



10~ 3 
= 0.1 . 



(24) 



In Figs. [2] and [3] respectively, we plot the numerical 
results for the evolution of the Hubble parameter and the 
equation of state. Also shown are zoomcd-in views of the 
evolution around the bounce point. One can see from 
Fig. [2] that the Hubble parameter H evolves smoothly 
through the bounce point with a dependence on cosmic 
time which is close to linear The maximal value of 



FIG. 3. Numerical plot of the evolution of the equation of 
state parameter w (vertical axis) as a function of cosmic time 
(horizontal axis). The insert shows the detailed evolution of w 
around the bounce time. The initial conditions were chosen 
to be: tj>i = —0.73 and 4>i = 1.5 x 1CT 4 . The background 
parameters were chosen as listed in Eq. ( |24[ |. Planck units 
are used. 



the Hubble parameter H , which we denote as the bounce 
scale Hb, is of the order O(10~ 4 M p ), and it is mainly 
determined by the value of the parameter Vq. 

From Fig. [3] one sees that the equation of state pa- 
rameter w of the scalar field is approximately equal to 
w = 5.67 which agrees w ell w ith what is obtained by in- 
serting q = 0.1 into Eq. |l8| [^] However, the Ekpyrotic 
contracting phase ends when H approaches the bounce 
scale Hb- We see that w then crosses the cosmological 
constant divide w = — 1, and this implies the violation of 
the Null Energy Condition in the bounce phase. After the 
bounce, the equation of state rapidly evolves back to be 
above the cosmological constant boundary and quickly 
approaches the value w = 1 which corresponds to the 
kinetic-driven phase of expansion. This fast-roll expand- 
ing phase is determined by the shape of the potential. 
During this period the contribution of the scalar field 4> 
will be diluted quickly relative to the contributions of 
regular matter and radiation. Thus, the universe in our 
model is able to connect smoothly to the usual thermal 
history of the Standard Big Bang model. 

We also plot the numerical results of the evolution of 
the background scalar field <p and its time derivative <f> 
as functions of cosmic time (see Fig. [4J. We see that 
the scalar 4> evolves monotonically from a large nega- 
tive value to a large positive one. Notice that far away 
from the bounce point, both the evolution of 4> and <\> are 
smooth, and <f> <C M p which implies that the higher order 



4 In the following section we make use of the linear approximation 
for H around the bounce point in our analytical study of the 
evolution of cosmological fluctuations. 



5 Since in the contracting phase the solution of is stable along 
with the background evolution, one can choose any arbitrary 
value of q to design a nonsingular bouncing model. 



7 



0.2 - 




0.12 - 
■ 0.08 - 
0.04 - 

0.00 * 

-100 -50 50 100 

t 

FIG. 4. Plot of the background scalar field (j> and its time 
derivative <j> as functions of cosmic time. The initial conditions 
and background parameters are the same as those chosen in 
Fig. [2] Planck units are used. 
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FIG. 5. Numerical plot of the ghost coefficient g as a function 
of cosmic time. The insert shows the detailed evolution of g 
close to the bounce point. During the interval when g > 1 
a ghost condensate will have formed. The initial conditions 
and background parameters are the same as those chosen in 
Fig. [2] Again, Planck units are used. 



operators in the Lagrangian ([I]) are highly suppressed by 
powers of M p . This explains why the approximate La- 
grangian ( 16 1 used in the semi-analytic analysis is valid. 
In the bounce phase, <j) suddenly obtains a dramatic en- 
hancement and then, after the bounce, falls back to a 
small value. Correspondingly, there is a sharp peak in the 
plot of <j) as a function of time. This is related to a short- 
term tachyonic instability about the bounce point which 
could lead to a "big rip" singularity as was found in the 
cosmology of the conformal Galileon [TTJtE]. However, 
since in the model the potential is bounded from below 
and the bounce phase only lasts a couple of Planck times, 
this instability does not have time to develop into some- 
thing which could destabilize the background solution. 

Finally, we calculate the evolution of the coefficients g 
and V (see Figs. [5] and [6]). From Fig. [5] We see that g <C 
1 far away from the bounce point, but becomes larger 
than one during the bounce phase which indicates that 
there exists a period during which a ghost condensate 
forms. However, the model is free of any ghosts in the 
far infrared since the sign of the coefficient V is always 
positive as shown in Fig. [6j 

To conclude this subsection, we have verified that our 
model leads to a nonsingular bounce which is free of 
instabilities at the level of homogeneous and isotropic 
cosmology. The stability of the model towards inhomo- 
geneities is investigated in the next section. 



F. On the Stability of the Bounce 

In recent work [46j . Steinhardt and Xue raised impor- 
tant concerns regarding the stability of modified Ekpy- 




FIG. 6. Numerical plot of the coefficient V as a function 
of cosmic time. The insert shows the detailed evolution of V 
near the bounce point. The initial conditions and background 
parameters are the same as those chosen in Fig. [2] Planck 
units are used. Note that the coefficient is positive throughout 
which implies the absence of ghost instabilities in the infrared 
limit. 



rotic bounces. They find that in the model [H] the 
anisotropy and the shear which are damped out during 
the phase of Ekpyrotic contraction shoot up to values 
larger than their initial values during the nonsingular 
bounce phase. In addition, they find that the adiabatic 
mode which is suppressed during the phase of Ekpyrotic 
contraction but which has a deep blue spectrum increases 
dramatically in amplitude during the bounce phase and 
dominates over the entropy mode which in the construc- 
tion of [U] has a scale-invariant spectrum. 



As we discuss in this subsection, our model appears 
to be free of these problems The model of [44] also 
uses a ghost condensate construction to achieve a non- 
singular bounce. However, the Null Energy Condition- 
violating ghost condensate appears at small values of X, 
namely for X < X C) where X c is some critical value. In 
the model of [H], the value of X at the beginning 
of the Ekpyrotic phase of contraction is larger than X c , 
and it increases rapidly during the Ekpyrotic contraction 
by a factor e 2N , where ./V is the number of e- foldings of 
modes which exit the Hubble radius during the Ekpyrotic 
phase. The relative contributions of the anisotropy and 
shear to the energy density decrease by this factor during 
the phase of Ekpyrotic contraction. In order to obtain a 
nonsingular bounce, X must decrease after the end of 
the Ekpyrotic phase by an even larger factor to reach the 
ghost condensate phase with X < X c < Xi. But this im- 
plies that the anisotropy and shear will increase by the 
corresponding factor and will come to dominate again. 
A corollary of this analysis is that the bouncing phase 
is long compared to the maximal value of the Hubble 
expansion rate. 

In our model, the ghost condensate is triggered not 
by X decreasing, but by </> increasing to some critical 
value. During the bouncing phase X can remain large. 
Hence, the bounce phase will be short on the time scale 
corresponding to the maximal value of the Hubble con- 
stant (this is verified in our numerical results), and the 
anisotropy and shear will not increase. Thus, it appears 
that our model is stable towards anisotropic stress and 
shear instabilities. 

Since the main idea of our model is to provide the 
adiabatic mode of curvature fluctuations with a scale- 
invariant spectrum, the concern of |46j regarding the 
spectrum of cosmological perturbations does not arise. In 
fact, the mode which we are following is the same mode 
as the one which becomes dominant after the bounce in 
[46) . However, in our case it has a scale- invariant spec- 
trum which is inherited from the matter-dominated phase 
of contraction which preceded the Ekpyrotic phase. We 
will return to this point at the end of the section on cos- 
mological perturbations. 



III. COSMOLOGICAL PERTURBATIONS 

We devote this section, to a study of the dynamics 
of linear cosmological perturbations in the model. At 
the linearized level, each Fourier mode of the fluctuating 
field evolves independently. It is useful to first consider 
a sketch of various relevant length scales in the nonsin- 
gular bouncing cosmology. The first length scale is the 
physical wavelength X p h — a/k of the fluctuation mode 



We are grateful to Paul Stcinhardt and BingKan Xue for detailed 
discussions on this point. 




FIG. 7. A sketch of the evolution of length scales in the non- 
singular bouncing universe. The vertical axis is the physical 
spatial coordinate x p h, and the horizontal axis is the cos- 
mic time t. The physical wavelength \ p h = a/k of the mode 
with comoving wavenumber k is depicted in green (the second 
curve from the bottom at the earliest time) ; the Hubble radius 
Ah = l-H] i s depicted in blue (the curve which diverges to 
infinity at the bounce point); and the red line (the constant 
bottommost line) denotes the Planck length A p = M" 1 . 



(labelled by comoving wavenumber k) which we wish to 
follow. This length must be compared with the Hubble 
radius H~ 1 {t). In order to allow for a causal generation 
mechanism of fluctuations, the wavelength must be sub- 
Hubble at very early times. The third length scale is the 
Planck length, the cutoff length below which our effective 
field theory breaks down. In Fig. [7] we present a sketch 
of the evolution of these length scales in the background 
cosmology. 

As is obvious from Fig. [7j the perturbation modes we 
are interested in originate inside the Hubble radius early 
in the contracting phase. They exit the Hubble radius, 
briefly re-enter at the bounce point, and finally enter the 
Hubble radius at late times in the period of Standard Big 
Bang evolution. We will separate the time evolution of 
the fluctuations into four phases: the initial era when 
the perturbations are set up (we will consider several 
possibilities in this section), the evolution in the Ekpy- 
rotic contracting phase, the dynamics during the ghost 
condensate-dominated bouncing phase, and the evolution 
in the post-bounce fast-roll expanding phase. We shall 
denote with subscripts, 'V', "b-"j "b" and "b+" the ini- 
tial moment when the perturbations are set up, the be- 
ginning of the bouncing phase, the bounce point, and the 
end of the bouncing phase, respectively. We will consider 
several choices for initial conditions for the cosmological 
perturbation: vacuum fluctuations during the Ekpyrotic 
phase, fluctuations formed in thermal equilibrium during 
the Ekpyrotic phase, and the fluctuations inherited from 
a period of matter contraction (the case we will be most 
interested in). 
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Since the equation of state is w > —1/3, the phys- 
ical wavelength X P h — a/k of a primordial fluctuation 
mode with fixed wavenumber k decreases less fast than 
the Hubble radius Xh = \tl\~ 1 . Thus, primordial per- 
turbations generated on sub-Hubble scales at the initial 
moment can escape into the super-Hubble regime during 
the contracting phase. At late times in the expanding 
phase, the Hubble radius catches up again and modes 
re-enter the Hubble radius. Modes which we are inter- 
ested in for cosmological observations today re-entered 
the Hubble radius at very late times (at times close to 
or after the time of equal matter and radiation). Thus, 
at the bounce point the wavelength is many orders of 
magnitude larger than the minimal value of the Hubble 
radius, which in turn is larger than the Planck length 
provided that the scale of the bounce is sub-Planckian 
(which is a necessary requirement for our effective field 
theory treatment to be justified). Thus, as in all non- 
singular bouncing cosmologies with bounce scale smaller 
than the Planck mass, fluctuations never enter the sub- 
Planckian region of ignorance (wavelengths smaller than 
the Planck length). This is a significant advantage com- 
pared to the situation in inflationary cosmology where 
the wavelengths of all fluctuation modes is smaller than 
the Planck length at the beginning of the phase of infla- 
tion if inflation lasted more than 70 e- foldings [57]. In 
our model, provided the energy scale of the universe at 
the bounce is lower than the Planck scale, then the phys- 
ical wavelength of a perturbation mode corresponding to 
the current Hubble radius is in the far infrared, as shown 
in Fig. [7j This lends strong support to the use of linear 
cosmological perturbation theory to study the evolution 
of the fluctuations. 



perturbations is to insert the ansatz for matter and met- 
ric into the action for matter and gravity and expand to 
quadratic order in the fluctuations about the cosmologi- 
cal background. Our detailed analysis is summarized in 
the first part of the Appendix. The resulting quadratic 
action for £ is: 



J drd 3 x Z ~ [C' : 



(25) 



where we have introduced conformal time r = J dt/a. 
The action is quadratic since we are working in linear per- 
turbation theory. The action contains no higher deriva- 
tive terms due to the special type of higher derivative 
Lagrangian under consideration. Thus, the only change 
compared to what is obtained in the case of a canonical 
scalar field coupled to Einstein gravity is the specific form 
of the speed of sound c s and of the function z which de- 
scribes the coupling between the background cosmology 
and the fluctuations. 

The parameter z is determined by the background met- 
ric and the background matter field. Its form is derived 
in the Appendix and given by Eq. ( 92 ) . The general 
expression of the sound speed parameter c s is provided 



by Eq. (93). In our model, these two parameters take on 
the form 

2M p 4 a 2 2 [M 2 (l - g) + 6 7 1^ + + ^} 



I 2 



(2M 2 H - 70 3 ) 2 
Af 2 (l-.g) + 4 7 i^ + /^ 2 -^g 
M 2 (l - g) + 6 7 i?0 + 3/3<j) 2 J- 5^ 



2 7 



(26) 



2 Mr, 



A. General Analysis 

Of the three families of metric fluctuations: scalar, 
vector and tensor modes, we will concentrate on scalar 
fluctuations, which couple at linear order in perturba- 
tion theory to matter. Vector perturbations and tensor 
fluctuations (gravitational waves) do not couple to scalar 
field matter at linear order. For matter given by a sin- 
gle scalar field there are - before fixing coordinates - four 
scalar metric degrees of freedom and one matter degree of 
freedom. However, two of these degrees of freedom can be 
eliminated by fixing coordinates. One degree of freedom 
is constrained by the Einstein constraint equation, and 
yet another vanishes since there is no anisotropic stress 
for scalar field matter. Hence, only one physical degree 
of freedom remains. We choose this variable to be the 
curvature fluctuation £ in uniform field gauge. 

To obtain the action for scalar cosmological perturba- 
tions, we make a gauge fixed ansatz for the metric and 
matter fluctuations, insert this ansatz into the Einstein 
equations and expand them to leading order (in the am- 
plitude of the fluctuations) about the cosmological back- 
ground. Another method of studying the dynamics of 



It is easy to check that for /? = 7 = g = the speed 
of sound becomes unity and the form of z reduces to 
the familiar one associated to a canonical scalar field in 
Einstein gravity [HI [33] . 

A common practice in the theory of cosmological per- 
turbations, is to introduce a convenient quantity related 
to C by, 



(27) 



which is a generalization of the well-known Mukhanov- 
Sasaki variable [SHI I n terms of this variable, the 



action ( 25 ) takes on the canonical form 



So = 



drd 3 a 



, (28) 



where a prime indicates a derivative with respect to con- 
formal time, and thus lends itself to the process of canoni- 
cal quantization. In Fourier space, the equation of motion 
for the Fourier mode is 



,2 7.2 



k'~ 



v k = 



(29) 
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B. Dynamics of Perturbations 

We discuss the solutions to the equation of motion for 
the cosmological fluctuations in each of the phases of the 
background evolution - the Ekpyrotic period, the bounce 
phase, and the kinetic-driven phase of expansion after 
the bounce. On the transition hypersurface between the 
phases we will continue the solutions making use of the 
matching conditions derived by Hwang- Vishniac [SD] , and 
by Deruelle-Mukhanov [61.. These conditions indicate if 
the induced metric on the matching surface and the ex- 
trinsic curvature are continuous. Note that we will be 
using matching conditions at the end of the Ekpyrotic 
phase of contraction, and at the end of the bounce phase. 
In both cases, the background dynamics also obey these 
matching conditions. Thus, the difficulties of matching at 
a singular surface between a contracting and an expand- 
ing phase which were discussed in |38] are not present in 
our case. This procedure is the same as the one used in 
the past (see e.g. [5J HI H] where the results of the 
approximate analytical calculations were also compared 
to direct numerical studies of the fluctuation equations) . 



1. Ekpyrotic contraction 

In the contracting phase, we have \<f>\ 1 and </> <C M p . 
Thus one obtains 5 — 0, and in addition the higher order 
operators in the Lagrangian are suppressed. In this case, 
we recover the same equations that apply for a normal 
canonical scalar field in Einstein gravity. In the current 
case the scalar field potential is an exponential function. 
As analyzed in the previous section, we obtain an Ekpy- 
rotic contracting phase and the background equation of 
state takes the value given in ( 18 ). In the approximation 
of equality in dT8b: 



Mia 



2„2 



, a oc (f B _ - t) !-« , 



(30) 



and Cg ~ 1. We have introduced the time moment f B _ 
when the scale factor would meet the big crunch singu- 
larity if there was no nonsingular bounce. If we were not 
interested in the bouncing phase, it would make sense 
to normalize the time axis such that f B _ = 0. In that 
case, we would find that the function g becomes unity 
slightly earlier, namely at a time ( 1 _ g ^ B (keeping in 
mind that T-Lb- is negative). This signals the beginning 
of the bounce phase. In our case, we choose the time axis 
such that r = is the midpoint of the bounce phase and 
t B— is the beginning of this phase. In this case 



TB- = T B - 



(1 



J? 

q)U B - 



(31) 



Therefore, the equation of motion for cosmological per- 
turbation in the contracting phase simplifies 



k 2 - 



q(2q-l) 



{l-qf{T-TB-, 



v k ~ 



(32) 



and the general analytical solution is 

v k( T ) = ci(k)y/r- f B -Ju c [k(r - f B - 



+C 2 (k) y/T-T B -Y Vc [fc(r - f B _ )] , (33) 



with 



l-'Sq 
2(1 -q) 



(34) 



and the subscript " c " indicates that we are discussing 
the solution in the contracting background as introduced 
in Sec. II D In Eq. (33), J and Y are two kinds of 
Bessel functions having indices v c . The coefficients C\ 
and C2 are only functions of comoving wave number k, 
and they are determined by the initial conditions of the 
cosmological perturbations which we will address in a 
subsequent subsection. For the moment, we keep the 
coefficients general. 

Using the small argument limiting form of the Bessel 
functions, the first mode in (33) is decreasing in time on 



super-Hubble scales while the second mode is constant. 
Since in the Ekpyrotic contracting phase the scale fac- 
tor decreases only very slowly, to first approximation, it 
remains true that the contribution to £ from the second 
mode dominates and it is therefore approximately con- 
stant in time. 

The initial power spectrum (see Subsection D) of 
the curvature fluctuations on super-Hubble scales can 
also be read off from the small argument expansion 
of the Bessel functions: The dominant mode scales as 
Y Vc {x) ~ x~ Va and hence the spectrum of the dominant 
mode of C, is 



P C (fc) ~ fc^"" |c 2 (fc)| 



2. Nonsingular bounce 



(35) 



When cj) evolves into the ghost condensate range, the 
kinetic term is no longer approximately canonical and the 
Null Energy Condition is violated. As we have discussed 
in the previous section, the universe will exit from the 
phase of Ekpyrotic contraction at some moment ts- , and 
the equation of state of the universe will cross w = — 1 
and decrease rapidly to negative infinity. Since in this 
period the 4> 2 term in the expression for the energy den- 
sity (JoJ) yields a negative contribution, it will eventually 
cancel all the other positive contributions to the energy 
density. This happens at a moment ts- We choose our 
time axis such that ts — 0. At this time, the Hub- 
ble parameter vanishes and the nonsingular cosmological 
bounce occurs. 

During the bounce phase, the deviation of the equation 
for fluctuations from the canonical one becomes impor- 
tant. When studying fluctuations in the bouncing phase, 
it has been shown in previous work [6HE] that it is a good 
approximation to model the evolution of the Hubble pa- 
rameter near the bounce as a linear function of cosmic 
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time 



H = Tt 



(36) 



where T is a positive constant which has dimensions of 
k 2 . This parametrization is valid in a class of fast bounce 
models, and the magnitude of T is usually set by the de- 
tailed microphysics of the bounce. In the specific example 
of o ur model with background parameter values from Eq. 
( [24b , T is of the order 0(1Q- 4 )M% (from Fig. [2}. 

IN ext we consider the determination of the sound speed 
square in the neighborhood of the bounce. From the 
expression (26), we find that, in addition to the Hubble 
parameter, one needs to know the form of (j) 2 around the 
bounce. We have evaluated this in our semi-analytical 
study of the background solution and the result is given 
in Eq. p2l. Combining Eqs. p2l, (p6l and po), we find 



that the sound speed parameter takes the approximate 
form 



-"sb 



V M4(g -1)2 



(37) 



b 



in Eq. (37) 



in the bouncing phase. The subscript 
indicates the bouncing phase. From this result, we see 
that the perturbations in the model are not completely 
stable around the bounce since it is possible for c 2 to 
be negative during the bouncing phase. For example, 



if we make use of the parameter choice (24), we obtain 
c 2 b ~ — |. A negative sound speed square leads to an 
exponential growth of the perturbation modes. We will, 
however, now show that this growth is not too large to 
loose perturbative control of the analysis. 

To calculate the growth of perturbations during the 
bounce phase, we return to the equation of motion ( 29 ) 
for the fluctuation modes. In that equation, the key 
quantit y is the parameter z whose explicit form is given 
by Eq. (26). In our model, <f> takes on its maximal value 



at the bounce point, as has been analyzed in Section II- 
D. The explicit expression is given in Eq. (22). When we 



choose the parameter Vo to be very small, we obtain the 
approximate value 



2M 2 (g - 1) 

3/3 



and thus we can simplify the expression for z 2 : 



(2M p 2 ff- 7 </> 3 ) 



:S\2 



(38) 



(39) 



in the bouncing phase. 

The dynamics of the z parameter during the bouncing 
phase depend mainly on the Hubble parameter H and 
the time derivative of the scalar <$>. We find that the 
evolution of <fi is approximated by 

(j) ~ <j) B e~^ , (40) 

where T is approximately one quarter of the duration 
of the bounce (as shown in Fig. pi. In the following 
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FIG. 8. Numerical result and analytical parameterization of 
<j> as a function of cosmic time around the bounce point. 



figure we compare this approximate expression with the 
numerical result for <fi. 

Inserting the above parameterizations for 



and (36) 

for H into the expression of z 2 , we can neglect the con- 
tribution from H and further simplify the result 



Therefore, we obtain 



z 
z 



2 3/?M p 4 

I T^- 



a 2 M 2 



with 



M 



T 



2 

J>2 



(2T 2 



6T 

J>2 



)t 2 



(41) 



(42) 



(43) 



Consequently, the perturbation equation ( 29 ) can be ap- 
proximated as follows: 



u 2 v k 







with 



-c 2 sb k 2 +a 2 B M 2 



(44) 



(45) 



in the bouncing phase. Here, as is the value of the scale 
factor at the bounce point which we set to unity. One 
of the solutions is exponentially growing, the other one 
decaying: 



v b k (r) c d 1 (k)e 



u(t-tb) 



d 2 e 



—ui(t—tb) 



(46) 



where d\ and d 2 are two coefficients which can be deter- 
mined by the matching conditions for fluctuations, and 
both of which inherit the spectrum of the dominant mode 
in the contracting phase. 

Note that for the infrared modes which we arc inter- 



ested in the first term in (45) is negligible. Hence, the 



growth rate of infrared modes during the bounce phase 
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is independent of k. The amplification factor T for the 
perturbation modes during the bounce phase is 



at the transition points between the phases, i.e. we must 



T 



rB+ 

oJ B- 



~ exp 



T 



2 2 



3TT 2 + T 2 T 



2rp4 



3T 4 JT 



2 

T 2 



B + 



B- 



(47) 



Inserting the values of the parameters T = 2.7 x 10~ 4 , 
T = 0.5, and t B+ = —t B _ = 1 (obtained from the numer- 
ical solution for the cosmological background) we obtain 



O(10 5 



(48) 



which is in approximate agreement with our numerical 
calculation as we demonstrate shortly. Note that this 
growth rate is the same for all infrared modes of interest 
for current cosmological observations. 



3. Fast-roll expansion 

After the bounce, the (p 2 term will decrease back to a 
value below the Planck scale and the scalar field <fi will 
cross the critical value <fi*+. Therefore, the phase of ghost 
condensation will cease at a moment t B + after the bounce 
and the universe will smoothly enter an expanding era. 
Afterwards, all higher order operators are suppressed by 
the Planck scale and thus the Lagrangian recovers its 
canonical form. Since the potential is very flat in the 
expanding era, the scalar field will, as discussed in the 
previous section, enter a fast-roll state with equation of 
state w ~ 1. In this period, the equation of motion for 
cosmological perturbations is given by 



1 



4(r - TB+) 



•i'k 



. 



and yields the following solution 



' ei{k)^/r~- 
+e 2 (fc) v / 7 



f B+ J [fc(r - f B+ )} 
- f B+ Y Q [k(T - T B +)] 



(49) 



(50) 



with f B + = t b+ 



2Hi 



The superscript "e" indicates 



that we are talking about the solution in the expanding 
phase. 

Modulo the square root factor, the first mode is con- 
stant whereas the second one is increasing logarithmically 
on super-Hubble scales. Hence, the second mode domi- 
nates at late times. In particular, we will be interested 
in the spectrum of that mode. 



C. Matching Conditions 

Having obtained general solutions to the perturbation 
equations in the various phases of the cosmological evolu- 
tion, we now study how the solutions are to be matched 



determine the coefficients appearing in the solutions ( 46 ) 
and (l50|) through the matching conditions. As mentioned 



at the beginning of this subsection, the matching con- 
ditions tell us that the induced 3-metric on the hyper- 
surface of the phase transition and its extrinsic curvature 
should be continuous across the matching surface. In an 
expanding universe, it has been shown [62U63] that in the 
absence of entropy fluctuations the curvature perturba- 
tion £ in constant field gauge is a conserved quantity on 
large scales (this result extends even beyond linear per- 
turbation theory [M]), and thus £ is the variable which 
is now generally used to describe the dynamics of cosmo- 
logical perturbations. However, in a contracting back- 
ground there are complications: £ has a growing mode, 
and care needs to be taken when applying the matching 
conditions. A second reason which forces us to reconsider 
the matching condition issue is that, when the universe 
experiences moments of violating the Null Energy Con- 
dition, i.e., w crossing w = — 1, the variable ( may not 
be well behaved since there are factors of (1 + w) which 
in the standard theory relate £ to other metric fluctua- 
tion variables. Therefore, we must consider the detailed 
dynamical evolution of the metric fluctuations in each 
phase and carefully study the matching conditions. 

In our model, we have chosen the uniform field gauge 
(see the first part of the Appendix). Since the field 4> 
is monotonically increasing in our scenario, this gauge 
is well-defined throughout, unlike what happens during 
inflationary reheating at a turning point of <f>(t). Our 
matching surfaces are the <fi — (/>*_ and <fi — 4>*+ sur- 
faces. The matching conditions for the background are 
satisfied on both matching surfaces (unlike what would 
happen if we were to try to match between a contracting 
branch and an expanding branch on a singular transition 
surface). Hence, the matching conditions of [60j [61] ap- 
ply, indicating that v and v' are continuous across the 
matching surfaces, i.e. at the moments t B _ and t B+ . 
Recall that both solutions of the perturbation equations 
on super-Hubble scales in both the contracting phase and 
kinetic-driven expanding phase are Bessel functions. We 
can use their asymptotic forms to simplify the calcula- 
tions. 

At the first transition surface, the transition between 
the contracting phase and the nonsingular bounce pe- 
riod, it is sufficient to consider the asymptotic form of 
the solution (33) in the contracting phase: 



Cl {k)k v - 
2 v «T Vc c 2 {k) 



(r-fB_) 3 ™c+ l 
{r-f B -Y 



(51) 



where T Vc is the Gamma function of z^ c -th order. To ob- 
tain this formula, we have made use of the expression for 
the index v c ( 34 ) and of the inverse of the relation of Eq. 



(181 to derive the explicit evolution of the fluctuation on 
super-Hubble scales as a function of conformal time. If 
the coefficients c\ and Ci are of the same order (whether 
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this is the case or not is determined by the initial condi- 
tions, as we shall discuss below) the second mode of v% 
(with coefficient function c 2 (k)) dominates at late stages 
of the contraction phase since the c\ term decreases faster 
as a function of time if w c is greater than 1. Q 

By requiring Vk and v' k to be continuous at the space- 
like surfaces of tb- and tb+, one can track how each 
Fourier mode evolves through the bouncing phase and 
derive the detailed expressions of all the coefficients ap- 
pearing in the above solutions of the perturbations equa- 
tions. We leave the detailed calculation for the second 
part of the Appendix, and here just write down the final 
expression of the perturbation which is 



T- 



n-i>c 



tt(t B - 



tb- 



TB+ 



TB. 



where 7e ~ 0.58 is the Euler-Masheroni constant. 

Let us comment on the time dependence and spec- 
trum of the resulting cosmological fluctuations, first fo- 
cusing on the amplitude. While the universe expands 
with an equation of state w = 1, the scale factor evolves 
as a ~ t2. This implies that the perturbation variable 
i>k evolves proportional to the scale factor and therefore 
the curvature perturbation £ will become conserved on 
super-Hubble scales after the bounce, as it must since 
there are no entropy fluctuations in the model. 

The spectrum shape is set by the spectrum of the co- 
efficient function 02(h) from the contracting phase. The 
late time power spectrum (see following subsection for 
a more detailed discussion) of the curvature fluctuation 
variable £, the power spectrum relevant for late time ob- 
servations, is proportional to the phase space factor P 
multiplied by the square of the absolute value of the mode 
function, i.e. 



P c (k) 



u3-2i/ c 



c 2 (fc)| : 



(53) 



Comparing with the spectrum of £ in the contracting 
phase ( 35 1 we see that the spectrum of ( has passed 



through the nonsingular bounce without change in its 
spectral index, but with boosted amplitude. This is the 
main result of this section. 



D. Initial conditions 

In this subsection, we investigate a group of initial con- 
ditions which set the momentum dependence of the co- 
efficient 02(h). We determine the conditions necessary to 



7 It is worth while to mention that the dominant mode is the ci 
term for a bouncing cosmology with w c less than 1, e.g. in the 
"Matter Bounce" _!/ cosmology. Since in general the cosmologi- 
cal perturbations after the bounce will have a spectrum which is 
determined mainly by the dominant mode of the fluctuations in 
the contracting phase, a different dominant mode will give rise 
to different results after applying the matching conditions. 



obtain a scale-invariant primordial power spectrum. Re- 
call the definition of the power spectrum of primordial 
curvature perturbations 

P P 

fjj , , ,11 ft, 



Pf 



2n 



\ Vk |2 



2tt 2 



p 



J k |2 



6ir 2 M? 



The spectral index is defined as 

dbxP, 



1 



C 



d\nk 



(54) 
(55) 

(56) 



Inserting the solution (52) into Eq. (54) and making 



use of the index ( 34 ) , we see that the primordial power 
spectrum will be scale- invariant (in the case w c > 1) 
when 



c 2 (h) ~ P c 



(57) 



is satisfieojf] 

We will first consider perturbations originating as vac- 
uum fluctuations in the phase of Ekpyrotic contraction, 
then perturbations which originate as thermal particle 
fluctuations, and finally the main case studied here, 
initial conditions inherited from vacuum perturbations 
which exit the Hubble radius in a matter-dominated 
phase of contraction which preceded the Ekpyrotic phase. 



1. Ekpyrotic vacuum fluctuations 

First, we consider cosmological perturbations which 
originate as quantum vacuum fluctuations on sub-Hubblc 
scales in the Ekpyrotic phase. Since the action for Vk is 
that of a harmonic oscillator, quantum vacuum initial 
conditions mean 



Vk(r) 



-ik(r— tb-) 



2k 



(58) 



To match this initial condition with the sub-Hubble per- 

(59) 



turbations in contracting phase (33), we get 

Ci ~ c 2 ~ fc° . 
Explicitly, we obtain 

c 2 ~ \ , (60) 

Inserting this result into the equation for the modes in 
the expanding phase, we obtain the final expression 



vl(r) 



Tj E a(r) 



2*k"°(T B - -f B - 



l a B + 



(61) 



8 In the case of w c < 1, e.g. as in the original "Matter Bounce" cos- 
mology, the ci mode dominates during the contracting phase. As 
a consequence, the matching conditions are modified and the con- 

_ _ 3 _ Guj c 

dition for scale invariance becomes c\ (k) ~ k v " 2 ~ fc i+3m c 
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Therefore, for vacuum initial conditions in the con- 
tracting phase, the power spectrum of primordial cosmo- 
logical perturbations becomes 



48tt 2 (1 + 3w c )t+^M2 \\hb 



(62) 



From this result, we can read that the largest amplitude 
of power spectrum is of order J- 2 ( ) 2 , the vacuum am- 
plitude of the perturbation mode which exited the Hub- 
ble radius at the time ig_ and for which k ~ \Wb- |j mul- 
tiplied by the boost factor J 72 . Modes with k > |%b-| 
will never exit the Hubble radius, they never undergo 
squeezing on super-Hubble scales, and hence remain as 
vacuum quantum fluctuations. Combining the above re- 



sult ( 62 ) with the definition of the spectral index from 
Eq. (56), we find that the spectral index for observable 



perturbation modes is 



7 + 9w c 
1 + 3w c 



(63) 



which is always blue tilted when w c is bigger than 
1 and converges to 3 when w c is large. This result 
agrees with previous studies of artificially smoothed 
out four-dimensional toy models of Ekpyrotic cosmology 
PS! [33 [65] 



2. Ekpyrotic thermal fluctuations 

In a bouncing cosmology, it is not manifest that pri- 
mordial perturbations have to arise as quantum vacuum 
fluctuations. An alternative choice for conditions for the 
primordial perturbations is thermal fluctuations of a gas 
of particles. If the fluctuations are generated by a ther- 
mal ensemble of point particles, e.g. relativistic particles 
with equation of state w r = 1/3, then it was shown that 
an Ekpyrotic contracting phase with a particular back- 
ground equation of state, namely w c ~ 7/3 is required to 
produce a nearly scale-invariant power spectrum [66] . In 
the following we will revisit this issue in the context of the 
explicit model of this paper. Our result is in agreement 
with the general analysis performed in Ref. |66j . More- 
over, we obtain the explicit relation between the spectral 
index and the model parameters. 

Consider that the universe is near thermal equilibrium 
at the beginning of the Ekpyrotic contracting phase. In 
that case, thermal fluctuations will be more important 
than vacuum perturbations. In a thermal system, the 
mean square mass fluctuation in a sphere of radius R 



9 Note that the fluctuations in Ekpyrotic models motivated from 
string theory may well be scale-invariant since the extra spatial 
dimensions }45| . entropy modes |40H44] and stringy effects (see 
e.g. [39]) play an important role. 



(and volume V) is given by the specific heat capacity 
C V (R) 



(SM 2 ) = C V (R)T 2 



(64) 



where T is the temperature of the thermal system. For 
a gas of normal particles, the heat capacity is extensive: 



C v — R 



dT 



R T 



(65) 



The mean square mass fluctuations in a sphere of radius 
R are given in terms of the k'th momentum mode (with 
k = aR^ 1 ) multiplied by the phase space factor fc 3 in the 
following way (see [33]) 



(66) 



where V = ^i? 3 . Inserting (64) and (65) into (66), one 
obtains 



Sp k ~ 0(l)fc-5i?-5T5 



(67) 



The density fluctuations determine the gravitational 
potential $ through the time-time component of the per- 
turbed Einstein equations (which is also the so-called rel- 
ativistic Poisson equation), and this in turn allows us to 
compute the Mukhanov-Sasaki variable Vk- The relations 
are 



k 



4nGSp k , 
a 2 H ( MpVk 



k 2 H 



(68) 
(69) 



while the length scale of the perturbation is smaller than 
the Hubble radius. As a consequence (making use of the 
fact that Vk oscillates in time with frequency k/a on sub- 
Hubble scales), the amplitude of the Mukhanov-Sasaki 
variable at Hubble crossing moment can be expressed as 



Vk(r H ) 



aHz k 8p k 
2iM^kH 



H=4 



ia k T k 
2e H M 2 k 2 



(70) 



where we have applied the Hubble radius crossing condi- 
tion k = aH. 

Since the temperature of the universe scales as T ~ a -1 
for normal radiation and since in the Ekpyrotic contract- 
ing phase a ^ r 1 ^, and since th ~ 1/fc is the Hubble 
crossing condition, we can derive the k dependence of the 
perturbation variable. The result is 



Vk{Tn{k)) ~ k 



which implies 

c 2 (fc) : 



HO 
2^r„„ 



k 2 v k (T H (k)) 



l-9-w 
k 2 (l + 3« 



(71) 



(72) 



Com paring Eq. ( 72 ) and the condition of scale invariancc 
(57), we conclude that only when w c ~ | the primor- 
dial power spectrum of curvature perturbations seeded by 
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FIG. 9. The spectral index of primordial power spectrum n 3 
(vertical axis) as a function of the parameter q (horizontal 
axis). The black dot corresponds to the case of exact scale 
invariance of the spectrum when q — 0.2. 



thermal fluctuations will be scale-invariant. This agrees 
with the result of |66j . 

Explicitly, the spectral index of the primordial power 
spectrum is determined by q as follows: 



(73) 



We plot the spectral index of the primordial power spec- 
trum seeded by thermal fluctuations as a function of the 
model parameter q in Fig. [9j From this figure, one can 
read off that when q is slightly less than 0.2, the primor- 
dial power spectrum would have a slight red tilt, a result 
which is favored by current CMB observations. 

To summarize, we have shown that, assuming ther- 
mal initial condition for cosmological perturbations in an 
Ekpyrotic phase of contraction, our model is able to gen- 
erate a nearly scale-invariant power spectrum by choosing 
a suitable parameter q. 



3. Fluctuations inherited from a phase of matter-dominated 
contraction 

If there is no matter and entropy generation during 
the bounce phase [^J then the universe must have orig- 
inated from a contracting phase which in the far past 
was dominated by matter and radiation. In particular, 
scales which we probe today in cosmological observa- 
tions are likely to have exited the Hubble radius during 



10 If the cosmological fluctuations remain in the linear regime, then 
one can show that no entropy is produced I67| . 



a matter-dominated contracting phase, and those which 
are probed by smaller scale observations (e.g. QSO ob- 
servations) would have been exited during a radiation- 
dominated contracting phase. Thus, it is logical to con- 
sider the spectrum of fluctuations in this context. Specif- 
ically, we must compute the spectrum of fluctuations for 
modes which have exited the Hubble radius in a pre- 
Ekpyrotic phase of contration. We will focus exclusively 
on adiabatic perturbations (given that we now have sev- 
eral matter fields, namely regular cold matter, radiation 
plus the scalar field </>, there is the possibility of entropy 
fluctuations) . 

To compute the spectrum of perturbations we make 
use of the gauge-invariant variable the curvature fluc- 
tuation in comoving coordinates, and the corresponding 
Mukhanov-Sasaki variable v = where z = \/2ea with 
e = —H/H 2 . The equation of motion for the Fourier 
mode Vk(r) in the context of standard Einstein gravity 
is 



v'l + (k 2 - —)v k = 



(74) 



where the prime denotes the derivative with respect to 
the comoving time as defined in the previous section. For 
a constant background equation of state w m = 0, one 
obtains 



With Ur„ = i- 



(75) 



We assume that the cosmological perturbations dur- 
ing the matter-dominated period of contraction originate 
from vacuum fluctuations, which implies that on sub- 
Hubble scales we must have 



1 



-i f T kdf 



2k 



(76) 



(for \kr\ 3> 1). This is cons iste nt with the vacuum initial 
condition provided in Eq. ( 58 1 when the last term is 
negligible. 

On super-Hubble scales we can use the small argument 
expansion of the Bessel function solution of Eq. ( 74 ) to 
find the asymptotic form 



T 2 



(77) 



(for \kr\ < 1). 

Matching the two asymptotic solutions (76) and (77) 
at the moment of Hubble crossing \kr\ ~ 1 yields the 
final form of the solution for vV l on super-Hubble scales 



(78) 



where f m — r m — 2 /rl m with r m defined as the end point 
of the period of matter contraction, and ri m is the co- 
moving Hubble parameter at that time. Thus, we have 
reproduced the well known result that the spectrum of 
curvature fluctuations originating from vacuum pertur- 
bations on scales which exit the Hubble radius during a 
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matter-dominated phase of contraction is scale-invariant 
O [35] (see also the detailed calculations in the Lee- Wick 
bounce model |68J). 

Since we have shown that the spectrum of curvature 
perturbations passes through the bounce phase with un- 
changed index, we conclude that on the scales consid- 
ered here, the final power spectrum of curvature fluctu- 
ations will be scale-invariant. Thus, we have shown that 
our current model provides a realization of the "matter 
bounce" scenario. 

To see this result in a more detail, we need to match the 



initial condition ( 78 1 with the asymptotical solution of Vk 



in the Ekpyrotic contracting phase (51) at the moment 
r m . This gives us the form of C2(k): 



c 2 (fc) 



Hr)fc"°-3(T m -f B -)"° 



(79) 



and thus we can derive the expression of the canoni- 
cal perturbation variable in the expanding phase of our 
model, the result being 



v k (T) ~ -T— — —air] 



(80) 



(52) 



where we have made use of Eq. 

As a consequence, the scale invariance of primor- 
dial fluctuations which exited the Hubble radius during 
the matter-dominated phase of contraction is preserved 
through the Ekpyrotic phase and the nonsingular bounc- 
ing phase, and the final power spectrum will be scale- 
invariant in the expanding phase. To find the amplitude 
of the spectrum, we apply the definition of P^ to obtain 



Pr 



T 1 



1%HI 
192tt 2 M2 



(81) 



where the parameter H m is the physical Hubble parame- 
ter at the end of matter-dominated period of contraction. 

Notice that if there is regular matter and radiation 
which dominated early in the contracting phase, the over- 
all spectrum of fluctuations will change its shape. It will 
be scale-invariant on very large scales (those which ex- 
ited the Hubble radius during matter contraction), it will 
then shift to being a vacuum spectrum for scales which 
exited in the radiation phase of contraction |69) , and will 



end up with the deep blue spectrum of ( 63 ) on scales 



which exit during the period of Ekpyrotic contraction. 



E. Numerical Analysis of Cosmological 
Perturbations 



Previously, we gave an analytic estimate of the am- 
plitude of primordial perturbations generated for vari- 
ous initial conditions in our model. The analysis in the 
bounce phase was rather non-trivial and the fluctuation 
modes undergo a period of exponential growth. In order 
to confirm the analytical approximations made and to 
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FIG. 10. Numerical plot of four groups of canonical perturba- 
tion modes Vk (vertical axis) as functions of cosmic time (hori- 
zontal axis). These modes are distinguished by the comoving 
wave number k, which are k — 10 (in red), k — 1CP 11 
(in orange), k — 10~ 7 (in green), and k — 10 -3 (in violet), 
respectively. The inner insert shows the detailed evolution 
of Vk during the bounce phase. The initial conditions of the 
background field and background parameters are the same as 
those chosen in Fig. 2. All numerical values are in Planck 
units M p . The initial conditions for perturbation modes were 
chosen as the vacuum initial condition given by Eq. (51). 



check that the model is indeed well-behaved when pass- 
ing through the nonsingular bounce point at the per- 
turbative level, we numerically studied the evolution of 
primordial perturbations throughout the bounce. Specif- 
ically, we considered the case of vacuum initial condition 
generated in the Ekpyrotic phase of contraction. 

First, we studied the time evolution of a set of canon- 
ical perturbation modes vj. as a function of cosmic time. 
The results are shown in Fig. 



10 The comoving wave 



numbers k chosen are (in Planck uni ts) 10 



-15 



10" 



10~ 7 , and 10 -3 . First, we see from Fig. 10 that the larger 
the comoving wave number is, the later the perturbation 
mode exits the Hubble radius and ceases oscillation (the 
purple mode is still oscillating). 

When the modes pass through the bounce point, their 
amplitude is amplified. For the long wavelength modes 
the amplification factor is clearly independent of k, and 
the amplitude of the amplification factor agrees well with 
our analytic estimate for T from Section III.B.2. From 
our graph, it is clear that the amplitude grows approx- 
imately exponentially with cosmic time. This is due to 
the fact that the parameter z undergoes a rapid change. 
The equation of motion for gravitational waves is similar 
to that of cosmological perturbations, except that the 
function z(t) is to be replaced by the scale factor a(t). 
Since the scale factor is approximately constant around 
the bounce point, the amplitude of gravitational waves 
does not jump at the bounce. Hence, the evolution of 
fluctuations about the bounce point produces a natural 
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FIG. 11. Numerical plot of the primordial power spectrum 
P( (vertical axis) as a function of comoving wave number k 
(horizontal axis). The initial conditions for the background 
field and background parameters are the same as those chosen 
in Fig. 2, and Planck units are used. The initial conditions 
for perturbation modes are the same as those chosen in Fig. 
|10| The spectral index is indicated in the plot: n s ~ 3. 



suppression of the tensor to scalar ratio of fluctuations by 
about a factor of F 2 which obviates the need for a mat- 
ter bounce curvaton. This growth of scalar modes which 
we find is closely related to the instability of nonsingular 
bounces discussed in the case of the New Ekpyrotic sce- 
nario in [46] . We have shown that in our bouncing model 
the instabilities of fluctuation modes remain under con- 
trol and do not destroy the predictions of the cosmolog- 
ical model. 

In Fig. [Tl] we show the scale dependence of the pri- 
mordial post-bounce power spectrum Pq as a function of 
the comoving wave number. The plot is on a log-log scale 
to accommodate the wide range of scales. From the fig- 
ure, we can read off that Pq is a power law function of 
k with the spectral index n s ~ 3 for perturbation modes 
directly generated from vacuum fluctuations. This result 
is exactly the same as what we obtained in Eq. (56) and 
the discussion performed there. 



IV. CONCLUSIONS AND DISCUSSION 

In this paper we have presented a single scalar field 
model which yields a nonsingular bouncing cosmology. 
It makes use of a negative exponential potential which 
yields an Ekpyrotic contracting phase. At high energy 
densities, the scalar field undergoes a phase transition 
to a ghost condensate. This leads to the violation of 
the Null Energy Condition which is required to obtain a 
nonsingular bounce in flat FRW models in General Rel- 
ativity. Following the bounce, a period of kinetic-driven 
expansion results. 

By adding regular matter and radiation to the model, 



we obtain a realization of the "matter bounce" scenario 
which is free of the anisotropy problem which plagues 
other realizations of this scenario. The cosmological sce- 
nario is therefore as follows: the universe begins in a 
contracting phase with cold matter dominating over rel- 
ativistic radiation and over the scalar field. Due to its 
equation of state w c > 1 the scalar field comes to dom- 
inate the universe and leads to an Ekpyrotic-type phase 
of contraction. This phase of contraction is free from the 
BKL instability since the energy density in anisotropies 
grows less fast than that in the scalar field. This phase 
is followed by the ghost condensate-driven bounce which 
in turn ends in a kinetic-driven expanding period. Even- 
tually, the usual radiation and matter come to dominate 
the energy density again, leading to the a Standard Big 
Bang expanding universe. 

We have performed a detailed study of the evolution of 
cosmological fluctuations in our model. We have shown 
that the spectrum retains its slope through the bounce. 
Thus, vacuum fluctuations which exit the Hubble radius 
during the matter-dominated phase of contraction ac- 
quire and maintain a scale-invariant spectrum. Thus, 
our model provides a realization of the "matter bounce" 
alternative to inflation which is free from the anisotropy 
problem which plagues previous realizations. We have 
also shown that in the absence of initial matter and ra- 
diation dominated phases of contraction, it is possible to 
obtain a scale- invariant spectrum of fluctuations from ini- 
tial thermal particle inhomogeneities provided the equa- 
tion of state during the contracting phase takes on a par- 
ticular value. 

We have found that the fluctuation modes undergo a 
period of exponential growth during the bounce phase. 
The growth factor T of the fluctuation mode during 
this phase, while being large in amplitude, is indepen- 
dent of k. Hence, the spectral shape passes through the 
bounce without change. The amplitude of scalar modes 
increases relative to that of tensor modes (for which 
T ~ 1). Hence, a small tensor to scalar ratio results. 
The increase of the amplitude of the fluctuations during 
the bounce phase has another implication: the value of 
the Hubble constant at the transition between matter- 
dominated contraction and Ekpyrotic contraction can be 
fairly low. If we compare the final amplitude (81) of 



the curvature power spectrum with what is required to 
match observations, then with J ~ 1 we would require 
a very large value of the Hubble expansion parameter at 
the beginning of the Ekpyrotic phase. In this context, 
the anisotropy problem might have reappeared: the ini- 
tial anisotropies cannot be larger than a critical value 
such that they begin to dominate at the beginning of the 
Ekpyrotic phase. 

In this paper we have focused on adiabatic fluctuations 
only. It would be interesting to study entropy modes in 
the model. 
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V. APPENDICES 

In the first part of this Appendix, we study the general 
treatment of linear cosmological perturbations within the 
KBG model. In the second part, we show how a Fourier 
mode of cosmological perturbations is transferred from 
the contracting phase to the bouncing phase and then to 
the expanding period through the matching condition. 



A. The derivation of the quadratic action of 
cosmological perturbations 

It is useful to study perturbation theory by making use 
of the ADM metric: 



ds 2 



N 2 dt 2 - h ij (dx i + N t dt)(dx j + N j dt) , (82) 



where N and N" 1 are the lapse function and shift vector, 
respectively. Making use of this metric ( [82] ) , one can 
decompose the original action of the model minimally 
coupled to Einstein gravity into time and space parts. 
To be explicit, the action can be written as 



S = / dtd 3 xVh, 



N 



M 2 (R 3 + KijK ij - K 2 ) 



+K{(j),X)+G{^X)U(j> 



(83) 



where R 3 is the Ricci scalar defined on the three- 
dimensional space and 



1 

2N 



(h tj - ViNj - VjNi) 



(84) 



is the extrinsic curvature. Since we want to investigate 
the cosmological perturbations which vary in space and 
time, one has to include the space-dependence in the ex- 
pression for the kinetic term of the scalar </> which is given 
by 



X 



1 

W 2 



{<(> - Nih ij dj 



1 



h ij di<pd j( f> . (85) 



We are interested only in the part of the action in- 
volving scalar metric and matter fluctuations. It is well 
known that, for a single scalar matter field minimally 
coupled to Einstein gravity, there exists only one scalar 



type of degree of freedom. We choose the uniform field 
gauge 







a e 



(86) 



so that the linear perturbations of the scalar field in the 
field Lagrangian can be eliminated. Further, the scalar 
components of the lapse function and shift vector can 
be determined through the Hamiltonian and Momentum 
constraints. Explicitly, the scalar contributions to the 
lapse function and shift vector take the form 



N =l + a , N t 



Mn 



(87) 



up to leading order. One can insert Eqs. (86) and (87) 



into the action ( 83 1 and then expand it to quadratic or- 



der. After a lengthy calculation, the form of the second 
order action is written as, 



S 2 = 



dtd 3 xa 3 



(2M 2 C - 2M 2 Ha + ^G, x a)^^ 



-3M 2 ( 2 - 2M 2 a 



dfC 



-M; 



3M 2 H 2 a 2 



6M 2 HaC-^ 6 G. x aC 



+ Y K ,xxa 2 + m^G^xa 2 



-(j) 2 {G^ + —G^ x <j>)a z 



-K*a 2 



-H^G^xxa 2 



which involves the perturbation variables £, a and a. 



Varying the action ( 88 ) with respect to a and a yields 



the following relations, 
2M 2 C 



2M£H 



d 2 a = 3M p a 2 ( 



3 G, X 

2M 3 d 2 C 

'^G.x 



(89) 



2MjU 



2M 3 a 2 ((-6M 2 H 2 



b 2 K, x + <j) 4 K :X x) 



(2M*H - <f> 3 G, x ) 2 
2M 3 a 2 ((12H l j> 3 G t x + 3H<j) 4 G ■ XX) 

{2M%H-&G tX ) 2 
2M p 3 a 2 C(20 2 G, + 4 G, X 0) 



(2M^H-^G, X ) 2 



(90) 



which are exactly the Hamiltonian and Momentum con- 
straints. 



Finally, by making use of the relations ( 89 ) , ( 90 ) and 



the quadratic action (88), we can obtain the action of 

n 

^ - -^t(^C) 21 



scalar perturbation at linear order, whose form is 

c 2 At) 



So = 



2 y ' 



a 2 (t) 



(91) 
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where we have introduced a function z which is com- 
pletely determined by the background evolution. Its ex- 
plicit form is given by 



2M^a 2 (j) 2 V 



(2M£H 



(92) 



with V being mentioned in Eq. (12), which is the coef- 



ficient of the second order derivative term of the back- 
ground equation of motion. From Eq. (92), one can 



immediately find that the positivity of the coefficient V 
directly determines the positivity of z 1 and thus can be 
used to judge whether there is a ghost mode or not. The 
square of the sound speed is given by 



G 



x< 



4Hcj)G,x 
+ (2G,x - 



A" 



2M2 



2G, 



't> 2 G, xx )4>] 



(93) 



B. Details of the matching condition calculation 



We first match the cosmological perturbation vi and v\ 
at the moment tb—, and then determine the coefficients 



d\ and d.2 as follows 



di 
d 2 - 



[l-KjjMfa--; f B -)] 



c 2 r„ 



<->(r B 



-■Kbjk v c(T B -—T B -) C 

T B-)[l_2iy c -2w(r B _ 



-TB-)] 



jk"c( TB _—f B -) 



(94) 
(95) 



in which the dominant contribution comes from the C2 
mode of the perturbation in the contracting phase. 

Similarly, we match the perturbations vi and v^. at 
the end of the bouncing phase tb+ and determine the 
coefficients e\ and e^. We make use of these coefficients 
and then can extract the dominant mode of w£ as follows, 



2 2 -^7r(r B _ - tb- 



xc 2 (fc)fc 



(96) 



Note that the boost factor T and the term inside the 
square brackets appearing in the numerator of this for- 
mula give the amplification factor which arises from the 
tachyonic instability when the perturbation mode evolves 
through the bouncing phase. This feature is quite in- 
teresting for phenomenological studies of cosmological 
perturbations in bouncing cosmologies. For example, in 
[70] a bounce curvaton scenario was proposed based on 
a tachyonic amplification of primordial isocurvature per- 
turbations which is very similar to the tachyonic insta- 
bility encountered here. The same process can help to 
improve the efficiency of preheating in a bouncing uni- 
verse [71] . Because of the large value of the amplification 
factor which we find in our current model, we do not need 
to invoke a bounce curvaton mechanism to suppress the 
tensor modes relative to the scalar one. 



[1] R. H. Brandenberger, "Introduction to Early Uni- 
verse Cosmology," PoS ICFI 2010. 001 (2010) 
[arXiv: 1103.2271] [astro-ph.CO]]. 

[2] J. -L. Lehners, "Ekpyrotic and Cyclic Cosmology," Phys. 
Rept. 465, 223 (2008) [arXiv:0806.T245| [astro-ph]]. 

[3] D. Wands, "Duality invariance of cosmological perturba- 



tion spectra ," Phys. Rev. D 60, 023507 (1999) [arXiv:gr 
qc/9809062_ 

[4] F. Finelli and R. Brandenberger, "On the generation of a 
scale-invariant spectrum of adiabatic fluctuations in cos- 
mological models with a contracting phase," Phys. Rev. 
D 65, 103522 (2002) [arXiv:hep-th/0112249| . 

[5] Y. F. Cai, T. Qiu, Y. S. Piao, M. Li and X. Zhang, 
"Bouncing Universe with Quintom Matter," JHEP 0710, 
071 (2007) [arXiv:0704.l"0 90 [gr-qc]]; 
Y. F. Cai, T. T. Qiu, J. Q. Xia and X. Zhang, "A 
Model Of Inflationary Cosmology Without Singularity," 
|arXiv:0808.0819| [astro-ph] . 

[6] Y. F. Cai, T. Qiu, R. Brandenberger, Y. S. Piao and 
X. Zhang, "On Perturbations of Quintom Bounce," 
JCAP 0803, 013 (2008) |ajXiv:0711.2l"87] [hep-th]]; 



Y. F. Cai and X. Zhang, "Evolution of Metric Perturba- 
tions in Quintom Bounce model," arXiv:0808.2551 [astro- 
ph]. 

[7] Y. -F. Cai, T. -t. Qiu, R. Brandenberger and X. - 
m. Zhang, "A Nonsingular Cosmology with a Scale- 
Invariant Spectrum of Cosmological Perturbations from 
Lee-Wick T heory," Phys. Rev. D 80, 023511 (2009) 
|arXiv:0810.4677] [hep-th]] . 
[8] C. Lin, R. H. Brandenberger and L. Levasseur Perreault, 
"A Matter Bounce By Means of Ghost Condensation," 
JCAP 1104, 019 (2011) jarXiv: 1007.2654] [hep-th]]. 
[9] S. Alexander, T. Biswas and R. H. Brandenberger, "On 
the Transfer of Adiabatic Fluctuations through a Nonsin- 
gular Cosmological Bounce," arXiv:0707.4679 [hep-th]. 

[10] T. Biswas, A. Mazumdar and W. Siegel, "Bouncing uni- 
verses in string-inspired gravity," JCAP 0603, 009 (2006) 
[arXiv:hep -th/0508194] . 

[11] T. Qiu, J. Evslin, Y. -F. Cai, M. Li and X. Zhang, 
"Bouncing Galileon Cosmologies," JCAP 1110, 036 
(2011) [arXiv:1108. 0593 [hep-th]]. 

[12] D. A. Easson, I. Sawicki and A. Vikman, "G-Bounce," 



20 



JCAP 1111, 021 (2011) [arXiv: 1109. 1 047 [hep-th]]. 

[13] P. Horava, "Quantum Gravity at a Lifshitz Point," Phys. 
Rev. D 79, 084008 (2009) [arXiv:0901.3775] [hep-th]]. 

[14] R. Brandenberger, "Matter Bounce in Horava- 
Lifshitz Cosmolo gy," Phys. Rev. D 80, 043516 (2009) 
[arXiv:0904.2835l [hep-th]] ; 

X. Gao, Y. Wang, W. Xue and R. Brandenberger, "Fluc- 
tuations in a Horava-Lifshitz Bouncing Cosmology," 
JCAP 1002, 020 (2010) |arXiv:0911.3196] [hep-th]]. 
[15] X. Gao, Y. Wang, R. Brandenberger and A. Riotto, 
"Cosmological Perturbations in Horava-Lifshitz Grav- 
ity," Phys. Rev. D 81, 083508 (2010) [arXiv:0905.382T1 
[hep-th]]; 

Y. -F. Cai and E. N. Saridakis, "Non-singular cosmology 

in a model of non-relativistic gravity," JCAP 0910, 020 

(2009) [arXiv:0906.l789] [hep-th]] . 
[16] R. Brandenberger, H. Firouzjahi and O. Saremi, "Cos- 
mological Perturbations on a Bouncing Brane," JCAP 

0711, 028 (2007) [arXiv:0707 4181] [hep-th]]. 
[17] D. A. Easson, R. Gregory, G. Tasinato and I. Zavala, 

"Cycling in the Throat," JHEP 0704, 026 (2007) [hepJ, 

th/0701252] [hep-th]] . 
[18] V. Sahni and A. Toporensky, "Cosmological Hysteresis 

and the Cyclic Universe," |arXiv7 l203.0395 [gr-qc]. 
[19] A. Kehagias and E. Kiritsis, "Mirage cosmology," JHEP 

9911, 022 (1999) [hep-th/9910T74] . 
[20] Y. Shtanov and V. Sahni, "Bouncing brane worlds," 

Phys. Lett. B 557, 1 (2003) |gr-qc/0208047] . 
[21] T. J. Battefeld and G. Geshnizjani, "Perturbations in 

a regular bouncing universe," Phys. Rev. D 73, 064013 

(2006) [hep-th/0503160] . 
[22] L. R. Abra mo and P. Pete r, "K-Bounce," JCAP 0709, 

001 (2007) |arXiv:0705.2893"l [astro-ph]]. 
[23] L. P. Chimento and M. I. Forte, "Anisotropic k-essence 

cosmologies, " Phys. Rev. D 73, 063502 (2006) |astro-| 

ph/05 10726] . 

[24] J. Martin and P. Peter, "Parametric amplification of met- 
ric fluctuations through a bouncing phase," Phys. Rev. 
D 68, 103517 (2003) [hep-th/0307077] . 

[25] D. M. Solomons, P. Dunsby and G. Ellis, "Bounce be- 
haviour in Kantowski-Sachs and Bianchi cosmologies," 
Class. Quant. Grav. 23, 6585 (2006) |gr-qc/01030"87] . 

[26] A. Anabalon and A. Cisterna, "Asymptotically (anti) de 
Sitter Black Holes and Wormholes with a Self Interact- 
ing Scalar Field in Four Dimensions," Phys. Rev. D 85, 
084035 (2012) [arXiv:1201.2"008] [hep-th]] . ' 

[27] J. Karouby and R. Brandenberger, "A Radiation Bounce 
from the Lee-Wick Construction?," Phys. Rev. D 82, 
063532 (2010) [arXiv:1004.4947| [hep-th]] ; 
J. Karouby, T. Qiu and R. Brandenberger, "On the Insta- 
bility o f the Lee- Wick Bounce," Phys. Rev. D 84, 043505 
(2011) [arXiv:1104.3193] [hep-th]] . 

[28] V. A. Belinsky, I. M. Khalatnikov and E. M. Lifshitz, 
"Oscillatory approach to a singular point in the relativis- 
tic cosmology," Adv. Phys. 19, 525 (1970). 

[29] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, 
"The ekpyrotic universe: Colliding branes and the origin 
of the hot big bang," Phys. Rev. D 64, 123522 (2001) 
|arXiv:hep-th/0103239| . 

[30] J. K. Erickson, D. H. Wesley, P. J. Steinhardt and 
N. Turok, "Kasner and mixmaster behavior in universes 
with e quation of state w > 1," Phys. Rev. D 69, 063514 
(2004) [hep-th/0312009] . 

[31] J. M. Bardeen, "Gauge Invariant Cosmological Pertur- 



bations," Phys. Rev. D 22, 1882 (1980). 
[32] V. F. Mukhanov, H. A. Feldman and R. H. Branden- 
berger, "Theory of cosmological perturbations. Part 1. 
Classical perturbations. Part 2. Quantum theory of per- 
turbations. Part 3. Extensions," Phys. Rept. 215, 203 
(1992). 

[33] R. H. Brandenberger, "Lectures on the theory of cosmo- 
logical perturbations," Lect. Notes Phys. 646, 127 (2004) 
[arXiv:hep-th/03060711 . 

[34] D. H. Lyth, "The Primordial curvature perturbation in 



the ekpyrotic universe," Phys. Lett. B 524, 1 (2002) hep 
ph/0106153 ; 

D. H. Lyth, "The Failure of cosmological perturbation 

theory in the new ekpyrotic scenario," Phys. Lett. B 526, 

173 (2002) [hep-ph/0110"007] . 
[35] R. Brandenberger and F. Finelli, "On the spectrum of 

fluctuations in an effective field theory of the Ekpyrotic 

universe," JHEP 0111, 056 (2001) |hep-th/0109004] . 
[36] J. -c. Hwang, "Cosmological structure problem in the 

ekpyrotic sce nario," Phys. Rev. D 65, 063514 (2002) 

|astro-ph/0109045"] . 
[37] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, 

"Density perturbations in the ekpyrotic scenario," Phys. 

Rev. D 66, 046005 (2002) hep-th/0109050 . 
[38] R. Durrer and F. Vernizzi, "Adiabatic perturbations in 

pre big bang models: Matching conditions and scale in- 

variance," P hys. Rev. D 66, 083503 (2002) [arXiv:hep-| 

ph/0203275] . 

[39] A. J. Tolley, N. Turok and P. J. Steinhardt, "Cosmologi- 
cal perturbations in a big crunch / big bang space-time," 
Phys. Rev. D 69, 106005 (2004) [hep-th/0306109]. 

[40] A. Notari and A. Riotto, "Isocurvature perturbations in 
the ekpyrotic universe." Nucl. Phys. B 644. 371 (2002) 
[hep-th/0205019| . 

[41] F. Fine lli, " Assisted contraction," Phys. Lett. B 545, 1 
(2002) [hep-th/ 0206~lT2] . 

[42] P. Creminelli and L. Senatore, "A Smooth bouncing cos- 
mology with scale invariant spectrum," JCAP 0711, 010 
(2007) [hep-th/0702T65] . 

[43] J. -L. Lehners, P. McFadden, N. Turok and P. J. Stein- 
hardt, "Generating ekpyrotic curvature perturbations 
before the big b ang," Phys. Rev. D 76, 103501 (2007) 
|hep-th/0702T53] [HEP-TH]] . 

[44] E. I. Buchbinder, J. Khoury and B. A. Ovrut, "New 
Ekpyrotic cosmo logy," Phys. Rev. D 76, 123503 (2007) 
[hep-th/0702154| . 

[45] T. J. Battefeld, S. P. Patil and R. H. Brandenberger, "On 
the transfer of metric fluctuations when extra dimensions 
bounce or stabil ize," Phys. Rev. D 73, 086002 (2006) 
|hep-th/0509043| . 

[46] B. Xue and P. J. Steinhardt, "Unstable growth of curva- 
ture perturbation in non-singular bouncing cosmologies," 
Phys. Rev. Lett. 105, 261301 (2010) arXiv:1007.2875 
[hep-th]]; 

B. Xue and P. J. Steinhardt, "Evolution of curvature and 
anisotropy near a nonsingular bounce," Phys. Rev. D 84, 
083520 (2011) [arXiv: 1 1 06 0416] [hep-th] ] . 

[47] N. Arkani-Hamed, H. -C. Cheng, M. A. Luty and 
S. Mukohyama, "Ghost condensation and a consistent in- 
frared modification of gravity," JHEP 0405, 074 (2004) 
[hep-th/0312099| . 

[48] A. Nicolis, R. Rattazzi and E. Trincherini, "The Galileon 
as a local modification of gravity," Phys. Rev. D 79, 
064036 (2009) [arXiv:0811.2197| [hep-th]]. 



21 



[49] D. B. Fairlie and J. Govaerts, "Universal field equations 

with reparametrization invariance," Phys. Lett. B 281, [61] 
49 (1992) |hep-th/9 202056 ; 

D. B. Fairlie and J. Govaerts, "Euler hierarchies and uni- 

versal equat ions," J. Math. Phys. 33, 3543 (1992) [hep^| [62] 

th/9204074] . 

[50] P. Creminelli, A. Nicolis and E. Trincherini, "Galilean 
Genesi s: An Alternativ e to inflation," JCAP 1011, 021 
(2010) [arXiv:1007.0"027| [hep-th]] . [63] 

[51] G. W. Horndeski, Int. J. Theor. Phys. 10, 363 (1974). 

[52] C. Deffayet, X. Gao, D. A. Steer and G. Zahariade, 

"From k-essence to generalised Galileons," Phys. Rev. [64] 
D 84, 064039 (2011) [arXiv:1103.3260l [hep-th]]. 

[53] C. Armendariz-Picon, V. F. Mukhanov and P. J. Stein- 

hardt, "A Dynamical solution to the problem of a small [65] 
cosmological constant and late time cosmic acceleration," 
Phys. Rev. Lett. 85, 4438 (2000) [astro-ph/0004134| ; 
C. Armendariz-Picon, V. F. Mukhanov and P. J. Stein- 
hardt, "Essentials of k essence," Phys. Rev. D 63, 103510 [66] 
(2001) astro-ph/0006373 . 

[54] C. Deffayet, O. Pujolas, I. Sawicki, A. Vikman, "Imper- 
fect Dark Energy from Kinetic Gravity Braiding," JCAP 
1010, 026 (2010). [arXiv:1008.0048] [hep-th]]. [67] 

[55] T. Kobayashi, M. Yamaguchi and J. 'i. Yokoyama, "G- 
inflation: Inflation driven by the Galileon field," Phys. 
Rev. Lett. 105, 231302 (2010) [arXiv: 1008 . 06031 [hep-th] ] . 

[56] J. -Q. Xia, Y. -F. Cai, T. -T. Qiu, G. -B. Zhao and 
X. Zhang, "Constraints on the Sound Speed of Dynami- 
cal Dark Energy," Int. J. Mod. Phys. D 17, 1229 (2008) 
|astro-ph/0703202] . [68] 

[57] R. H. Brandenberger, "Inflationary cosmology: Progress 
and problems," |hep-ph/991Q410} 

R. H. Brandenberger and J. Martin, "The Robustness of 

inflation to changes in superPlanck scale physics," Mod. 

Phys. Lett. A 16, 999 (2001) |astro-ph/0005432| ; [69] 

J. Martin and R. H. Brandenberger, "The TransPlanck- 

ian problem of inflationary cosmology," Phys. Rev. D 63, 

123501 (2001) [hep-th/0005209] . 
[58] M. Sasaki, "Large Scale Quantum Fluctuations in the In- [70] 

flationary Universe," Prog. Theor. Phys. 76, 1036 (1986). 
[59] V. F. Mukhanov, "Quantum Theory of Gauge Invariant 

Cosmological Perturbations," Sov. Phys. JETP 67, 1297 [71] 

(1988) [Zh. Eksp. Teor. Fiz. 94N7, 1 (1988)]. 
[60] J. c. Hwang and E. T. Vishniac, "Gauge-invariant joining 

conditions for cosmological perturbations," Astrophys. J. 



382, 363 (1991). 

N. Deruelle and V. F. Mukhanov, "On matching condi- 
tions for cosmological perturbations," Phys. Rev. D 52, 
5549 (1995) |arXiv:gr-qc/ 9503050 . 

J. M. Bardeen, P. J. Steinhardt and M. S. Turner, "Spon- 
taneous Creation Of Almost Scale - Free Density Pertur- 
bations In An Inflationary Universe," Phys. Rev. D 28, 
679 (1983). 

R. H. Brandenberger and R. Kahn, "Cosmological Per- 
turbations In Inflationary Universe Models," Phys. Rev. 
D 29, 2172 (1984). 

D. Langlois and F. Vernizzi, "Evolution of non-linear cos- 
mological perturbations," Phys. Rev. Lett. 95, 091303 
(2005) [astro-ph/0503416| . 

S. Tsujikawa, R. Brandenberger and F. Finelli, "On the 
construction of nonsingular pre - big bang and ekpyrotic 
cosmologies and the resulting density perturbations," 
Phys. Rev. D 66, 083513 (2002) [hep-th/0207228] . 
Y. -F. Cai, W. Xue, R. Brandenberger and X. - 
m. Zhang, "Thermal Fluctuations and Bouncing Cos- 
mologies," JCAP 0906, 037 (2009) [arXiv:0903.4938| 
[hep-th]]. 

R. H. Brandenberger, T. Prokopec and V. F. Mukhanov, 
"The Entropy Of The Gravitational Fie ld," Phys. Rev. 
D 48, 2443 (1993) |arXiv:gr-qc/9208009| ; 
R. H. Brandenberger, V. F. Mukhanov and T. Prokopec, 
"Entropy of a classical stochastic field and cosmolog- 
ical perturbations," Phy s. Rev. Lett. 69, 3606 (1992) 
[arXiv:astro^ph7 9206005 . 

Y. -F. Cai, T. -t. Qiu, R. Brandenberger and X. - 
m. Zhang, "A Nonsingular Cosmology with a Scale- 
Invariant Spectrum of Cosmological Perturbations from 
Lee-Wick Theory," Phys. Rev. D 80, 023511 (2009) 
[arXiv:0810.4"677| [hep-th]]. 

H. Li, J. Q. Xia, R. Brandenberger and X. Zhang, 
"Constraints on Models with a Break in the Primor- 
dial Power Spect rum," Phys. Lett. B 690, 451 (2010) 
|arXiv:0903.3725l [astro-ph.CO]]. 

Y. -F. Cai, R. Brandenberger and X. Zhang, "The Mat- 
ter Bounce Curvaton Scenario," JCAP 1103, 003 (2011) 



arXiv: 110 1.0822 [hep-th]]. 
Y. -F. Cai, R. Brandenberger and X. Zhang, "Preheat- 
ing a bouncing universe," Phys. Lett. B 703, 25 (2011) 
[arXiv:1105.4"286l [hep-th]] . 



